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ABSTRACT 

We assess the influence of massive black hole (MBH) ejections from galaxy centres 
due to gravitational radiation recoil, along the cosmic merger history of the MBH 
population. We discuss the 'danger' of recoil for MBHs as a function of different MBH 
spin-orbit configurations and of the host halo cosmic bias, and on how that reflects on 
the occupation fraction of MBHs. We assess ejection probabilities for mergers occurring 
in a gas-poor environment, in which the MBH binary coalescence is driven by stellar 
dynamical processes and the spin-orbit configuration is expected to be isotropically 
distributed. We contrast this case with the 'aligned' case. The latter is the more 
realistic situation for gas-rich, i.e., 'wet,' mergers, which are expected for high-redshift 
galaxies. We find that if all haloes at z > 5-7 host a MBH, the probability of the Milky 
Way (or similar size galaxy) to host a MBH today is less than 50%, unless MBHs form 
continuously in galaxies. The occupation fraction of MBHs, intimately related to halo 
bias and MBH formation efficiency, plays a crucial role in increasing the retention 
fraction. Small haloes, with shallow potential wells and low escape velocities, have 
a high ejection probability, but the MBH merger rate is very low along their galaxy 
formation merger hierarchy: MBH formation processes are likely inefficient in such 
shallow potential wells. Recoils can decrease the overall frequency of MBHs in small 
galaxies to ~ 60%, while they have little effect on the frequency of MBHs in large 
galaxies (at most a 20% effect). 
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1 INTRODUCTION 



■ It is well established that galaxies contain massive black 
5-H ' holes (MBHs) in their nucle i (e.g.. iKormendv fc Richstond 
5^ 1 1 19951 : iRichstone et alj Il998t ). Observations of the host el- 
liptical galaxy or spiral bulge show a relation between 
the mass of the MBHs and the galacti c spheroid lumi- 
nosity and the stellar velocity dispersion jMagorrian et al.l 
1 19981 : iFerrarese fc MerrittfeOOd ; iGebhardt et al.ll2000l) . Such 
relationships, as wel l as their small in strinsic scatter 
(|Tremaine et al.ll2002l : iGiiltekin et all 120091 1. hint at a con- 
nection between the formation of the MBH and the forma- 
tion of the spheroid. 

In the context of the favored cold dark matter cosmol- 
ogy, hierarchical assembly of galaxies and protogalaxies with 
MBHs at their centers naturally lead s to the formation of 
MBH binaries (|Begelman et al.l Il980t ). If the MBH binary 
can achieve a small enough separation through stellar dy- 
namical hardening or through viscous gas dynamical drag, 
gravitational radiation from the binary can become signif- 
icant enough to cause the system to merge. Gravitational 



waves carry away energy from the binary system, driving the 
system to merger. If the black holes have unequal masses or 
misaligned spins, anisotropic gravitational radiation will im- 
part a net momentum flux on the binary, causing the center 
of mass to recoil. 

Until recently, the magnitude of the recoil ve- 
locity was uncertain. For non-spinning black holes 
the maximum recoil ha s now been calc ulated to be 
Wrecoil.max ~ 200 km s -1 |Baker et al.l 120061 ) . and a sim- 
ilar range is expected for black holes with low spins, 
or with spins ( ant i-) aligned with the binary orbital an- 
gular momentum. Several studies have also found con- 
sistent results for the maximum recoil of spinning black 
holes t> r ecoii.max = 1000-4000 kms~* dBaker et al j 120071; 
Gonzalez et al.l 2007 : Herrmann et al.l 2007 : Koppitz et al 



20071 : ICampanelli et alj 120071 : | Schnittman fc Buonannol 
20071) . Such velocities are significant because they will eject 
the binary from the galaxy since most galaxies have escape 
speeds < 1000 km s -1 . 

IVolonteri et alj (|2008l ) studied the effect of recoil on 
the MBH occupation fraction in nearby galaxies. They as- 



© 0000 RAS 



2 Volonteri et al. 



sumed that the relative orientation between the orbital an- 
gular momentum of the binary and the spins of the two 
MBHs were isotropically distributed. This configuration can 
resul t in high recoil velo c ities, up to thousands kms -1 . How- 
ever, iBogdanovic et all (|2007l ) proposed that MBHs orbit- 
ing in gaseous circumnuclear discs, such as th ose expected 
in ad vanced stages of gas rich galaxy mergers (|Maver et al.l 
2007), align their spins with the orbital angular momen- 
tum of the binary. This configuration leads to small recoils 
for the MBH remnant. The accreting gas exerts gravito- 
magnetic torques that suffice to align the spins of both the 
MB Hs with the angula r momentum of the large-scale gas 
flow. iDotti et al.l |2009h have quantified the efficiency of this 
alignment process through the analysis of high resolution TV- 
body Smoothed Particle Hydrodynamic s (SPH) simulations . 
They apply the algorithm presented in IPerego et all (|2009l ^ 
to evolve masses, magnitudes and orientation of the spins of 
two MBHs. 

We extend the investigation of IVolonteri et al.1 (|2008l ) 
considering different degrees of alignment between the spins 
and the angular momentum of the MBHs. We will either 
use isotropic spins and angular momentum configurations 
or the 'quas i - align ed' distributions of spin orientations from 
IDotti et al.1 (|2009l V Those two cases bracket all the possi- 
ble values of recoil velocities. We apply these distributions 
to Monte-Carlo realizations of galaxy merger trees, and we 
study the effect of gravitational recoil on the occupation 
fraction of MBHs in galaxies at different redshifts. 



2 RECOIL VELOCITIES 
2.1 Fitting formulae 

The recoil velocity may be broken down into components 
arising solely from mass asymmetry, v m , which is perpen- 
dicular to the orbital angular momentum vector L pa i r , and 
components arising from spin asymmetry, v± and , which 
are perpendicular and parallel to L pa i r , respectively: 



^'recoil 



^Jv^ + vj_+ 2v m vj_ cos(£) + u 2 , 



(1) 



where £ is the angle between v m and v ± in the orbital plane 
We as s umed £ = 145° , as s uggest ed by 
' 2009h . ICampanelli et al.l (|2007l ) and 



Lousto & 


z Zlochowei 


Lousto & 


i Zlochower 



20091 , fit CL) ropose the following fitting formulas for the 
recoil components: 

v m = ArfsJT^ 4^(1 + B77), (2) 
v± = Hrt^l + qy'ial-qal), (3) 

t)|| = Kr/ 2 (1 + cos(0 — 60) Ui" — ga^j , (4) 

where q — M2/M1 ^ 1 is the mass ratio of the black holes 
and Tj = q/(l + q) 2 is the symmetric mass ratio. The compo- 
nents of the spins of the two MBHs are broken into projec- 
tions parallel and perpendicular to L pa j r , a ± — as'm(6) and 
a" = acos(#). In the isotropic case, cos(#i) and cos(#2) are 
distributed uniformly between —1 and 1. We compare the 
results from the isotropic case with those obtained assuming 
quasi-aligned spin-orbit configurations. In this case we as- 
sumed 6\ and 82 were distributed according to IDotti et al.l 
|2009l . see Section [2.2[) . Q is the angle between (a^ — qo.\) 



and the separation vector at coalescence, and Oo depends 
on the initial separation between the holes. We assume a 
uniform distribution of — Go between and 2ir. Note that 
for qa\± = a-2i_ (including a\± — ei2± = 0), vu vanishes 
and there is no recoil out of the orbital plane. In the ab- 
sence of spins, the recoil is v m , which is maximized for for 
q = (3 — \/5)/2 ~ . 38. T he best fit parameters found by 
lLousto fc Zlochowerl (|2009l . which we refer to as fit CL) are 
A = 1.2 X 10 4 kms"\ B = -0.93, H = 6900kms-\ and 
K = 6.0 x 10 4 kms~ 1 

iBaker et al.l l)2008l . referred to as fit B) found an alter- 
native form for v\\ that scales as 77 s : 



Vu = Krj 3 (l + q) 1 ( af cos(<J>i) — qa^ cos($2j 



where $1 and $2 are the differences between two angles 
that depend on the system of reference. We as sume <3?i = 
$2 u niformly distributed between and 2-zr (see Dotti et al.l 
l2009h . The best-fit parameters from lBaker et al.1 (|2008f ) are 
A = 1.35 x 10 4 kms _1 , B = -1.48, H = 7540 kms -1 , and 

K = 2.4 x 10 5 km3" ' . 

iHerrmann et al.l (|2007l . fit H in the following) parame- 
terize the recoil velocity in terms of 9h, the angle between 



(6) 



where the / are the spin angular momenta of the holes. Tak- 
ing Lp a i r to be in the z direction, they find the Cartesian 
components of the recoil velocity to be: 



V x = C H x cos{6n), 
V y = C H y cos(en), 
T4 = Co-K*sin(0H), 



(7) 



where C Q = Eg 2 (Afi + M 2 )' 2 (l + g) -4 with the best-fit 
parameters H x = 2.1 X 10 3 , H y = 7.3 x 10 3 , and K z = 
2.1 x 10 4 . 

In Figure [T] we plot recoil velocities from fits H, CL, 
and B as a function of q for a sample of isotropically dis- 
tributed MBH spins and for aligned configurations. This fig- 
ure demonstrates the agreement of the three formulae for dif- 
ferent spin-orbit configurations. For each value of q between 

0. 1 and 1, we produce 500 000 realizations of the two MBH 
spins with magnitudes uniformly distributed between and 

1. For the isotropic case we also assign isotropically dis- 
tributed spin directions. We do not show results for q < 0.1 
because the recoil velocity drops below the typical escape 
velocity from galaxies. 

Fit B decreases faster as q decreases, with respect to 
fits CL and H. This is because the component of the recoil 
velocity parallel to L pa i r scales as rf in fit B, and as rj 2 in 
fits CL and H. The average values of « rcco ii obtained with 
the different prescriptions for isotropically distributed spins 
are consistent within la (i.e. within a factor of ~ 2). As- 
suming perfect alignment between the two spins and L pa i r , 
the results from fits CL, H, and B are consistent within a 
factor of 2.5 for any q. For the aligned case with q fts 0.1, 
the difference between fits CL and B disappears because, for 
perfect alignment, a = a" and «y = 0, resulting in a perfect 
agreement between the two fits. 
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Figure 1. Recoil velocity of the MBH remnant as a function of q, 
for MBH spins isotropically distributed (top) or for aligned con- 
figurations (bottom). The upper, middle, and lower panel refers 
to fit H, fit CL, and fit B, respectively. The thick lines show the 
average values of u recoi i (jJ reC oil), and tJ r ccoil ± a. 



2.2 Quasi-aligned configurations 

As mentioned in the Introduction, gravo-magnetic torques 
exerted by accreting flows onto the MBHs tend to align a\ 
and a 2 to Lpair- The distributions of relative angles between 
the two MBH spins, and between each spin and ipair af- 
ter the formation of a MBH binary have been computed in 
iDotti et all (|2009l , see Figure 2 therein). Their results stem 
from the analysis of high resolution Af-body/SPH simula- 
tions. In the suite of runs discussed in that paper, the Au- 
thors varied the thermodynamical properties of the gas. In 



particular they assumed two different polytropic equations 
of state, with polytropic index 7 = 5/3 ('hot' runs) and 7/5 
('cold' runs), respectively. The hot case corresponds to an 
adiabatic monoatomic gas, as if radiative cooling were com- 
pletely suppressed during the merger. This case m imics gas 
radiatively heated by an AGN (jMaver et al.ll2007l ). The cold 
case, instead, has been shown to provide a good approxima- 
tion to a gas of sola r metallicity heated by a starburst (e.g., 
ISpaans fc ' Silkll200Ch. 

In each run, iDotti et al] (120091 ) found a significant de- 
gree of alignment already present before the formation of 
a binary. However, their runs always start with an equal- 
mass MBH pair. Because the dynamical evolution of the 
two MBHs is faster than the Salpeter time, the masses of 
the two MBHs do not change significantly, and at the end of 
the runs a nearly equal-mass binary forms. Here we discuss 
the depend e nce o f spins/orbit configurations discussed by 
IDotti et all (|2009h on the masses of the MBHs, by compar- 
ing the dynamical timescale (rd yn ) of the binary formation 
with the alignment timescale (r a ii g n). Before we consider the 
ratio of these two timescales, we look at each in slightly more 
detail in t he context o f our s imulations, but for a full dis- 
cussion see IDotti et al.1 (|2009l ). 

Dynamical friction is the process that drives the dynam- 
ical evolution of the pair. Since the longer timescale is what 
determines the time until coalescence, Td yn is the dynamical 
friction time-scale of the less massive MBH, which depends 
on a number of factors: 



Tdy 



(8) 



where K is a function of th e properties of the circumnuclear 
disc. In IDotti et~aH (|2009T ) runs, Mbh = 4 x 10 6 M Q , and 
rdyn = 4.5 — 7.5 Myr, depending on the initial orbital param- 
eters of the BHs and on the effective equation of state used 
to evolve the gas thermodynamics. H ere we are consider ing 
the circumnuclear disc parameters in IDotti et al.1 (|2009r ) a.s 
fiducial. Considering different values for these is out of the 
scope of this paper, but it is likely to change only the details 
and not our qualitative results. 

The alignment timescale also depends on the BH spin, 
a, and may be expressed in terms of the Eddington fraction 
(see equation 43 in Perego et al. 2009): 



Tali, 



n 5 / 7 M~ 2 / 35 f~ 32 / 35 
1 BH JEdd 



(0) 



For Mbh = 4 x 10 6 Mq in the simple case of a fixed orbital 
plane of the outer disc, i.e., coherent ac cretion, the v alue o f 
Taiign ~ 10 yr. For the simulations in IDotti et al.1 l|2009h . 
however, the orbital plane of the accreted particles is never 
constant, and thus the alignment timescale will be longer: 
Taiign ~ 1-4 Myr, depending on the initial orbital parame- 
ters of the BHs and the effective equation of state used. 

Since both black hole spins must achieve alignment with 
the disc angular momentum for the spins to be aligned with 
each other, the alignment timescale we are interested in 
is the longer of the two. For most cases, this will be the 
timescale for the less massive binary. Before the formation 
of a binary, the more massive of the two MBHs always ac- 
cretes more mass, and has the shortest alignment timescale 
(|Dotti et al.1 I2009T ) . Thus, we may use Mbh = M 2 in the 
expression for T a iig n . 

The situation after the formation of a hard binary, how- 
ever, is slightly different. In this case, the torques exerted by 
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the binary onto the gas carve out a low-density region in the 
middle of the disc, and the secondary MBH, being closer to 
the outer and denser gas region, can grow faster than the 
primary (jHavasaki et al.ll2007l ; ICuadra et alj|2009h . further 
aligning its spin with L pa i r . For our purposes, we conser- 
vatively assume that after forming a binary, the spin-orbit 
configuration at the time of binary formation does not sig- 
nificantly evolve. 

The ratio of these two quantities is obtained by dividing 
Eq [9] by Eq [8] and shows the dependence on the accretion 
behaviour of the two MBHs: 



Tali, 



Jr dy 



Mo 



2/35 r-32/35 
JEdd 



(10) 



For a fixed value of a and rewriting in terms of /Edd 
MM 2 _1 , this expression becomes 



7"align/7dyn 



M. 



-32/35^13/7 



(11) 



If they are accreting at the Bondi rate, M2 oc M 2 , and 
consequently Taiign/Tdyn °c M.\ ' 3B . In this case the ratio is 
almost independent of the mass of the secondary, and the 
assumption that the two MBHs are almost aligned is always 
correct. 

If, instead, the seco ndary is accr e ting a t the Eddington 
limit, Mi oc Mi. From iDotti et all (|2009l ) we know that 
for M 2 « 4 X 10 6 Mq, r a i ign /r dyn « 0.30 Using this as 
our normalization, we get r a ii gn /rdy n ~ O.lMj 3 / 35 , where 
A/2,6 is the mass of the secondary in units of 10 6 Mq. This 
scaling implies that a binary can form before the two MBHs 
align their spins to L pa i r only for very massive secondaries 
(Mi > 2 x 10 7 Mq) and rapidly accreting pairs. 

To estimate the relevance of such mergers, we con- 
sider the prevalence of the secondary MBH having mass 
M2 > 2 x 10 7 M Q . Assuming the scaling of MBH masses with 
the velo city dispersion of th e host determined for nearby 
galaxies (|Giiltekin et alibOQSl ). the host of a MBH with mass 
M2 > 2 x 10 7 Mq has a > 130 km s _1 , and an escape velocity 
> 1000 km s -1 . The recoil velocity has a significant proba- 
bility of being > 1000 km s -1 only for nearly equal-mass 
mergers (q > 0.3). Such merg ers are extremely rare for the 
heav iest MBHs (see section 4; IVolonterll2007l : ISesana et al.l 
2007). Furthermore, as mentioned above, further accretion 
can increase the degree of alignment of the binary system. 
As a consequence, in this study we may neglect the depen- 
dence of the spin-orbit configurations on the MBH masses 
and adopt the same distributions of relative angles between 
the two MBH spins for all MBH mergers in our sample. 
Note, also, that the isotropic distribution case we test be- 
low provides an upper limit to the ejection probability at all 
masses. 



3 EJECTIONS AND HALO BIAS 

In models of structure formation based on gravitational in- 
stabilities in Gaussian primordial fluctuations, the number 
density and bias properties of a halo can be expressed as a 
function of its rms fluctuation of the linear density field at 



1 We took the cold, prograde case as our fiducial case, but the 
results are not very different for the other cases. 



redshift z, a(Mh,z), and on the threshold density for col- 
lapse of a homogeneous spherical perturbation at redshift z, 
8 c (z). 

In the lPress fc Schechterl (|l974T ) formalism the comov- 
ing number density of haloes of mass between M and 
M + dM can be expressed as: 



dn 
dM 



2 prn 
7T U 



-d(lno-) 
dM 



(12) 



where v c = S c (z)/a(M, z) is the number of standard devi- 
ations which the critical collapse overdensity represents on 
mass scale M. At any redshift we can identify the charac- 
teristic mass (i.e., v c = 1), and its multiples. The higher v c , 
the more massive and rarer the halo, and the higher its bias 
and clustering strength. Although this formalism is derived 
from a Press-Schechter analysis, it agrees fairly well with the 
resu lts of TV-body simulat i ons. 

IMartini fc Weinbergl (|200ll ) suggested that clustering 
analysis of quasar samples can be deconvolved to yield the 
typical mass of haloe s hosti ng quasars, and their typical v c . 
IMartini fc Weinberg (l200lf) and sub s equent investigations 
(|Shen et al.1 120071 : iMvers et all 120071 : IPorciani et"alTl2004r i 
have found that high redshift quasars are highly biased ob- 
jects with respect to the underlying matt er, and that their 
increases with z. IShen et ail (|2007l) and IMvers et al.l (|2007f) 
find that the bias increases from v c ~ 3 at z = 2, to v a ~ 3.5 
at z = 3, and v c ~ 5.5 at z = 5. 

We evaluate the ejection probabilities, using the tech- 
niques described in section 2 with fit CL, for haloes repre- 
senting peaks of the density fluctuations u c =l, 2, 3, 4, 5, 6 
as a function of reds hift in a concordance ACDM cosmology 
(|Spergel et al.l |2007| 1 . For every halo mass we estimate the 
ejection probability by comparing the recoil velocity (v rcco n) 
to the escape velocity from the dark matter halo potential 
well (truncat ed at the virial radiu s) . We model the halo po- 
tential with a lNavarro et al.l (|l997l ) de nsity profile, where th e 
halo properties evolve as suggested bv lBullock et al.l (|200lh . 

The ejection probabilities thus defined are a lower 
limit to the probability that the recoil voids a galaxy of 
its central MBH, as if the recoil velocity is lower than 
the escape velocity but higher than the velocity disper- 
sion of a halo, the timescale for the recoiled MBH to re- 
turn to the center under the effect of dynamical friction 
is likely longer than the Hubble time ( Madau fc QuataertJ 
2004 iGualandris fc Merrittll2008l: IVolonteri fc Madaull200l 



Devecchi et al.l 120091 : iGuedes et ail 120091 . and references 
therein) . 

With this caveat in mind we can estimate lower limits to 
the probability that MBH binaries are ejected or displaced 
due to the gravitational recoil in ^ c -peaks haloes (Figure [2]). 

The redshift at which the probability drops to 50% in- 
creases with increasing v c for all binary mass ratios and 
spins. The probability of ejection for MBH binaries in haloes 
which host very high redshift quasars {v c ~ 5.5 at z = 5) 
drops to at z > 19 if MBHs are non-spinning, while for 
spinning MBHs the ejection probability is significant down 
to z > 13 in the absence of any alignment process. Lower v c 
peaks remain dangerous for merging binaries for longer. A 
v c ~ 3 peak has an ejection probability of 100% until 2 = 13. 
If all v c ~ 3 haloes host a MBH, and all experience a major 
merger before z = 13, then no MBHs are left for further 
evolution. In this scenario, MBHs must form at z < 11, in 
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Figure 2. Ejection probability as a function of redshift for 
different u c = 8 c (z)/cr(M, z) peaks of the density fluctuations 
field. Line styles indicate the assumed spin of the merging black 
holes: Solid lines: a = 0.9; Dashed lines: a = 0.6; Dotted lines: 
a = 0. Line colors indicate assumed distribution of spin orienta- 
tions: Green: isotropically distributed; Red: orientation distribu- 
tion from hot simulations by Dotti et al. (2009); Blue: orientation 
distribution from cold simulations by Dotti et al. (2009). All prob- 
abilities assume q = 0.1. 



order to evolve all the way to z = 3 and host the observable 
quasars. 

However, the j/ c -peak formalism does not allow a clear 
evaluation of the merger history of haloes. A halo which rep- 
resents a i/ c -peak at a given redshift z, would be incorporated 
into a halo representing a lower peak at z — Az. A MBH in a 
given halo would then belong to different ^ c -peaks in its life- 
time. In order to evaluate the number of MBH mergers that 
a halo experienced it is necessary to use techniques that can 
trace the whole merger history of a given halo as a function 
of mass and time. 



4 INDIVIDUAL HALO HISTORIES 

We now turn to evaluate the merger histories of different 
galaxy haloes. We adopt a statistical approach, based on 
merger histories extracted from full ACD M merger trees. 
This w ork has a more general approach than lVolonteri et al.l 
l|2008l ). and addresse s in an astroph ysical context th e anal - 
ysis by ISchnittmanl (|2007l ). While IVolonteri et~aH (|2008h 
looked at the role of ejections in the merger histories for 
galaxies that were in clusters but not in the 'main trunk,' 
here we examine merger histories with a focus on the main 
halo of the merger tree. We evaluate the average ejection 
probability as a function of cosmic epoch, and consider 
how it depends on the environmental bias. We consider 
here the halo merger histories leading to the formation of 
haloes with masses M =4x 10 13 M Q , M = 2 x 10 12 M©, 
Mo = 2 x 10 11 Mq at 2 = 0. We average our results over 
20 realizations of the same mass, to account for cosmic vari- 
ance. Our technique a nd cosmological frame work is similar 
to the one described in lVolonteri et all (|2003l ). We track the 
dynamical evolution of MBHs ab-initio and follow their as- 
sembly down to z = 0. 



Several theoretical arguments indicate that MBH for- 
mation occu rs at very high redshift, and probably in biase d 
haloes (e.g.. iMadau fe Rees|[200ll ; IVolonteri fe Reesll2006t ). 
Additional parameters, such as the angular momentum of 
the gas, its ability to cool and its metal enrichment, are 
likely to set the exact efficiency of MBH formation and the 
redshift range whe n the mechanism ope rates (for a thor- 
ough discussion see IVolonteri et al]|2008l ). We consider the 
effect of varying the initial conditions, by assuming frequent 
or rare MBH formation process. As a reference model, we 
assume that MBH formation is effective in all haloes with 
v c > 3 at z > 20 (corresponding to masses > 10 6 Mq). This 
model is based on a scenario where MBHs form as remnants 
of the first generation of metal-free stars (Population III). 
The main features of the hierarchical assembly of MBHs left 
over by th e first stars in a ACDM cosmology have been dis - 
cuss ed by IVolonteri etail (120031), IVolonteri fe Reesl (|2006t ) 
and IVolonteri fe Natarai an (2009). A model similar to the 
one presented in this paper has been shown to reproduce ob- 
servational constraints on MBH evolution (luminosity func- 
tion of quasars and Soltan's argument, Mb h-j relationship 
at z = 0, mass density in MBHs at z = 0; IVolonteri et al.l 
120081 ). We also analysed a very different case, where when- 
ever a halo grows above a given mass threshold (> 10 10 M ) 
and it does not already contain a MBH in its centre, then 
it forms one (in case the previous MBH had been displaced 
or ejected, a new MBH materializes), regardless of redshift. 
This model is not based on a specific physical model, but it 
is meant to be used for comparison with existing numerical 
simulations. MBH formation mechanisms define when and 
how often haloes are populated with MBHs. They provide 
the initial occupation fraction. 

Additionally we have to follow the dynamical evolution 
of haloes and embedded MBHs all the way to z = in order 
to determine the effect of recoils on the occupation fraction 
of MBHs. Since the magnitude of the recoil depends on the 
mass ratio of the merging MBHs, we have to model the mass- 
growth of MBHs and the merger efficiency. We base our 
model of MBH growth on the empirical correlation found be- 
tween MBH masses and the properties of their hosts, and on 
the suggestion that these correlations are established during 
galaxy mergers that fuel MBH accretion and form bulges. 
We therefore assume that after every merger between two 
galaxies with a mass ratio larger than 1 : 10, their MBHs 
climb to the same relation with the velocity dis persion of 
the halo, as it is seen today (|Giiltekin et al.l 20091 ) : 



M = 1.3 x 10 M M ( 







200 kms- 



(13) 



We link the correlation between the black hole mass 
and the central stellar velocity dispersion of the host with 
the empirical correlation between the central stellar veloc- 
ity dispersion and the asymptotic circular velocity (V c ) of 
galaxies (Ferrarese 2002; see also Pizzella et al. 2005; Baes 
et al. 2003). 



200 kms" 



304 kms- 



(14) 



The latter is a measure of the total mass of the dark 
matter halo of the host galaxies. A halo of mass Mh collaps- 
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ing at redshift z has a circular velocity 



aker Campanelli Herrmann 



V c = 142 km s" 



1 


M h 


1/3 




A c " 


1/6 




_10 12 M M0 _ 




O 2 


18tt 2 _ 





(l + z) 1/2 

(15) 

where A c is the over-density at virialization relative to the 
critical density. For a WMAP5 cosmology we adopt here the 
fitting formula (Bryan & Norman 1998) A c = 18tt 2 + 82<f - 
39d 2 , where d = — 1 is evaluated at the collapse redshift, 

so that = n ro (i + z ) 3 /(n m (i + z f + n A + n fc (i + «) 2 ). 

Therefore, if a dark matter halo of mass Mh at red- 
shift z hosts a MBH, we can derive the MBH mass via V c 
and a (see al so IVolonteri et~all2003l ; iRhook fc Wvithell2005l ; 
ICroton|[2009l ). The relationship between black hole and dark 
matter halo mass is: 



M = 1.3xlO 8 M 



M h 


5/3 




A c 1 


9.41O 13 M _ 




Oz 
. ^ 'm 


18tt 2 _ 



-i 5/6 



(l+*) 5/2 , 



in goo d agreement with the recent result bv lBandara et all 
l|2009h . 

We further assume that MBHs merge within the merger 
timescale of their host haloes (t merge ), which is a likely as- 
sumption fra_^ffil^biiraries formed after gas rich galaxy 
mergers |Dotti et al.l 120071 . an d references therei n). Wc 
adopt the relations suggested bv lTaffoni et al.l |2005t ) for the 
orbital decay of merging haloes. We treat as MBH mergers 
at a given z those MBH binaries that are expected to merge 
during that specific timestep. If two haloes start interacting 
at Zi n , corresponding to a Hubble time tn (zin) with a merger 
timescale t mer ge, then we consider these MBHs merged at Zfi n 



corresponding to tff(z fln ) = t H (zin) + t v 



Although we 



follow the dynamical evolution of each MBH along cosmic 
time, we note that dynamical friction appears to be efficient 
for mergers with mass ratio of the progenitors larger than 
1 : 10. Small satellites suffer severe mass losses by the tidal 
perturbations induced by the gravitational field of the pri- 
mary halo. This progressive mass loss increases the decay 
time that can be of order the Hubble time when the mass 
ratio is smaller than 1 : 10. 

Finally, the spin of MBHs is fixed to be a = 0.9, in order 
to obtain upper limits to the recoil consequences. 

This set of models, though simple, provides a wide range 
of histories that can be considered to bracket some extreme 
behaviours. We find that the environment of the MBH pop- 
ulation plays an important role. Along the merger history 
of massive galaxies, the fraction of ejected MBHs decreases 
more rapidly with decreasing redshift, dropping below 50% 
by z ~ 7 for M = 4x 10 13 M Q . The 50% threshold is reached 
at later times in the 'trees' of less massive haloes: z ~ 5 for 
Mo = 2 X 10 12 Mq and z ~ 2 for M = 2 x 10 11 M Q . 

It is particularly instructive to examine the ejection his- 
tory of the tree's main halo, which is the galaxy we would 
see today (Figure [3]). The main halo is usually among the 
most massive haloes in the tree at any time, implying the 
largest escape velocities. The fraction of ejected MBHs is 
therefore lower. The central MBH in a large galaxy has a 
very small probability of being lost after z — 5. Since most of 
the growth of MBHs happens between z — 3 and z = 1, cor- 
responding to the peak in quasar activity, we argue that, if 
accretion is responsible for creating the correlations between 
MBH s and their hosts (|Silk fc Reesil 19981 : iFabianll 19991 : iKind 
l2003h . then MBHs hosted in large galaxies will sit close to 




5 10 15 5 10 15 5 10 15 
redshift redshift redshift 



Figure 3. Probability of MBH ejections as a function of redshift 
along the merger history of the main halo in each merger tree. 
Each row corresponds to different halo masses Mo at z = as in- 
dicated in the left panels. Each column as indicated corresponds 
to the different fitting function used to calculate recoil speed. The 
green, vertically hatched regions are for an isotropic distribution 
of spin orientations; the red, horizontally hatched regions are for 
the 'hot' simulations by Dotti et al. (2009); and the blue, diago- 
nally hatched regions are for the 'cold' simulations by Dotti et al. 
(2009). In every panel the MBH spins are fixed at a = 0.9. The 
trends and results are in strong qualitative agreement with the 
results from the simpler simulations presented in Fig. [2] There is 
also general agreement among the different fitting functions. 



the expected correlation. The central MBH in a small galaxy 
(e.g., Mo = 2 x 10 11 Mq) has instead a large probability 
(~ 20%) of being ejected all the way to today. If such a 
low-redshift ejection happens, it is not immediately clear if 
the galaxy can re-acquire a MBH and if the latter can grow 
to the hypothetical mass that the correlation with the host 
would suggest. 

The alternative case of MBH formation (haloes with 
> 10 10 M o ) gives quantitatively similar results for large 
galaxies, e.g., Mo — 4 x 10 13 Mq, and qualitatively simi- 
lar results for small galaxies. Black holes form much later in 
this model (as the emergence of more massive galaxies must 
wait until z ~ 5 — 7), and this pushes the ejection fraction to 
first rise steeply, as more and more MBHs are available to 
merge, and then decrease as galaxies grow in mass deepening 
the potential wells. 

It is most interesting that the cosmic ejection rate is 
very similar for fit B compared to fits CL and H, notwith- 
standing the different r\ dependence (rj i vs rf 2 ). The reason 
is simple. From Figure [T] it is evident that the mass-ratio 
dependence picks up at q < 0.1 (where, incidentally, re- 
coil velocities become much smaller than escape velocities 
from gala xies). However, mergers with q <C 1 are ra ther un- 
common (IVolonteri et al ] 120031 : Is csana et al .1120051 ). due to 
dynamical effects. During a galactic merger, it is dynami- 
cal friction that drags in the satellite, along with its central 
MBH, towards the center of the more massive progenitor. 
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When the orbital decay is efficient, the satellite hole moves 
towards the center of the more massive progenitor, leading 
to the formation of a bound MBH binary. The efficiency of 
dynamical friction decreases with the mass ratio of the merg- 
ing galaxies: only nearly equal mass galaxy mergers ('major 
mergers', mass ratio larger than ~ 1 : 10) lead to efficient 
MBH binary formation within timescales shorter than the 
Hubble time (|Taffoni et al ] |2003h . These effects must be con- 
volved with the mass-ratio probability distribution. As the 
mass function of haloes (and galaxies) is steep, the proba- 
bility of halo mergers decreases with increasing mass ratio. 
That is, dynamically efficient major mergers are rare, and 
minor mergers are common but inefficient at forming MBH 
binaries. 

We can now derive how strong or mild recoils influ- 
ence the frequency of MBHs in galaxies. We analyse here 
the v c > 3 case only, as by construction the alternative 
case (mass threshold > 1O 1O M0) has an occupation frac- 
tion of unity for all haloes above the threshold, as if a MBH 
is lost, a new immediately materialises. . We first derive 
a control case, in which we ignore recoils altogether. This 
control case is shown in the bottom panel of Figure [4] The 
MBH frequency is calculated above a given minimum halo 



(> 10 10 M ? ; > 10 11 M w ; > 10 



12 Ma ) As discussed 
by iMenou et alj (120011 ) and IVolonteri et al.l (|2003h . the fre- 
quency of MBHs in haloes above a fixed mass threshold ini- 
tially decreases with cosmic time as lower mass haloes lack- 
ing MBHs become more massive than the assumed thresh- 
old. Eventually, the occupation fraction starts to increase as 
the total number of individual haloes drops. 

Given that different fitting formulae give very similar 
results, we show here the results that we obtain using fit CL. 
The worst case scenario, the case with the largest number 
of ejections, is the isotropic case, which leads to the largest 
recoils (see Figure Q]). The top panel of Figure [4] indicates 
that recoils can decrease the overall frequency of MBHs in 
small galaxies to ~ 60%, while they have little effect on the 
frequency of MBHs in large galaxies (at most a 20% effect). 
The middle panel depicts the cold aligned case. We consider 
the latter the more realistic situation for gas-rich mergers, 
which are expected for high-redshift galaxies, at least. It is 
reassuring to notice that recoils are not dangerous for most 
galaxies. 

This is the result of a combination of effects: on the 
one hand MBHs hosted in small haloes, with shallow po- 
tential wells and low escape velocities, have a high ejection 
probability (Figure [3]), on the other hand the MBH merger 
rate is very low along their galaxy formation merger hierar- 
chy: MBH formation processes are inefficient in such shal- 
low potential wells, and the anti-hierarchical nature of the 
galaxy assembly implies that not much action happens in 
low-bias systems at high-redshift. This is exemplified in the 
right panels of Figure [4] in which we show the frequency of 
MBH close binaries for the same mass thresholds we used 
to calculate the MBH frequency (note the different y-axis 
scales). Naively, one would expect the frequency of double 
MBHs to scale as the square of the MBH frequency, but the 
frequency of close binaries is suppressed with respect to the 
frequency of pairs because of the requirement of efficient or- 
bital decay. Our simulations show that gravitational recoil 
is not expected to be efficient at ejecting MBHs with mass 
~ 10 6 Mq . This mass range is where the Laser Interfero- 




2 4 6 
redshift 



Figure 4. Left: frequency of MBHs in galaxies as a func- 
tion of redshift for haloes above three different mass thresholds: 
M > 10 10 M Q (gray, vertically hatched) M > 10 11 M (magenta, 
diagonally hatched) M > 10 12 Mq (cyan, horizontally hatched) 
for isotropic spin orientations (top) and spin orientations from 
the 'cold' simulations of Dotti et al. (2009) (middle). The bottom 
panel shows the black hole occupation fraction when the recoil is 
arbitrarily set to zero at each MBH merger. Right: binary MBH 
frequency for the same models and mass thresholds. 



metric Space Antenna (LISA) gravitational wave observa- 
tory will be mo st sensitiv e to extreme mass ratio inspiral 
(EMRI) events (|Gairll2009l ). Thus the event rate for EMRIs 
is not likely to be affected by recoils. 

Summarising, the larger ejection probability that MBHs 
hosted in low-bias haloes have is counteracted by their lower 
merger probability, thus leading to an overall small change 
in the MBH frequency as a function of the recoil strength. 



5 DISCUSSION 

Ejection of MBHs may explain the unusual case of 
NGC 3115. NGC 3115 is an SO gala xy with MBH measured 
to ha ve mass M = 1 X 10 9 Mq (|Kormendv fc Richstonej 
yet the surface-brightness profile of the central bulge 
has a small break radius r& ~ 2 pc with a stee p inner slope 
slope consistent with a power-law 7 = 0.52 (|Lauer et al.l 
120071 ). This presents a contradiction of theoretical expecta- 
tions for a galaxy of this size (Mh ~ 10 13 Mq). 

A galaxy as massive as NGC 3115 is expected to have 
merged frequently enough to have had a binary b lack hole at 
some point in its history (|Volonteri et a l. 2003). According 
to th e standard pi c ture o f core scouring (jBegelman et al.l 



.Faber et al.1 Il997l ; iMilosavlievic 
Lauer et alj|2007l ). the binary black hole would eject stars 



198(j 



fc Merrittl |2003| : 



on elongated orbits as orbital energy is transferred from the 
binary to kinetic energy in the stars. This process shrinks 
the binary separation and flattens the stellar number den- 
sity profile so that the surface brightens profile appears as 
a core at the center. The galaxy's core can be described by 
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its deficit in stellar mass, i.e., the mass in stars ejected from 
what was previously a power-law profile. The most recent 
numerical simulations on black hole mergers find that for 
nearly equal masses, the mass deficit sh ould scale with the 
total mass of the binary |Merrittl [200^ 1 ; but for mass ra- 
tios far from unity, th e mass deficit shou ld scale with the 
mass of the secondary jSesana et al-lliooih . Observationally, 
the mass deficit is f ound to scale linearly with e xisting Mbh 
(|Lauer et alj 120071 ; iKormendv fc Benderl [2009h ■ Regardless 
of its magnitude, the mass deficit is expected to persist af- 
ter the binary coalesces because the remnant black hole acts 
as a 'guardian' tha t prevents dynamical refilling of the core 
(|Lauer et alj|2007h . 

The observations of NGC 3115, however, reveal that 
this is not what happened. There is an existing MBH at 
its center that would have acted as the core's guardian, 
but the surface brightness profiles show no core. There are 
several possible explanations. First, core scouring may be 
greatly diminished if the final major merge r had a signifi- 
cant amount of gas (|Kormendv et aLll2009h . The gas drag 
speeds up the process of MBH merger before the MBH bi- 
nary can eject many stars. Even if core scouring occurred, 
a subsequent gas-rich merger could re-fill the core by de- 
positing gas that would form stars. Here, we suggest another 
possibility. The final MBH binar y is ejected by gravitational 
recoil (or three-body en counters, iGiiltekin et al. |2004| . 120061 ; 
iHofhnan fc Loebll2007h . but the MBH is replaced in a sub- 
sequent merger. Taking into consideration the typical mass 
ratios of merging MBHs as a function of galaxy bias and 
cosmic time, the probability that a MBH in NGC 3115 has 
been ejected at redshift z < 5 is 10% in the case of a 'dry 
merger', but less than 1% in the case of a 'wet merger'. This 
is because at late cosmic times the binary mass ratio distri- 
bution becoming shallow, with q <C 1 becoming less probable 
(|Volonterill2007f ). In general, in order to have a probability 
larger than 10% that a MBH recoils with a velocity larger 
than the escape velocity of the progenitor of NGC 3115 at 
z ~ 5 the mass ratio of the merging MBHs must be q > 0.3 
for the isotropic and hot cases (assuming spins a = O.9.). For 
the cold case, even for q = 1 and a = 0.9, the probability 
is 0.8%. In the ejection scenario NGC 3115 should have had 
at least one major merger that involved two MBH-hosting 
galaxies where substantial star formation or AGN feedback 
kept the gas pressurized (corresponding to our hot simula- 
tions) . 



6 SUMMARY 

We assessed the effects of spin alignment on the strength of 
the gravitational recoil along the cosmic build-up of galaxies 
and MBHs. We determined ejection probabilities for merg- 
ers in gas-poor galaxies, where the MBH binary coalescence 
is driven by stellar dynamical processes, and the spin-orbit 
configuration is expected to be isotropically distributed. We 
contrast this case of gas-rich mergers, where we expect MBH 
spins to align with the orbital angular momentum. This is 
because in gas-rich environments MBHs accrete gas, which 
exerts gravito-magnetic torques that align the spins of both 
the MBHs with the angular momentum of the large-scale 
gas flow. We find that for aligned configurations the ejection 
probability is strongly suppressed (by at least a factor of 2). 



Along the merger history of a large elliptical, the ejection 
probability becomes negligible at z < 5, while small galax- 
ies have ejection probabilities of order 20% even today. How- 
ever, the MBH merger rate is very low along their merger 
hierarchy of small galaxies: MBH formation processes are 
likely inefficient in such shallow potential wells. The occu- 
pation fraction of MBHs, intimately related to halo bias and 
MBH formation efficiency, therefore plays a crucial role in 
increasing the retention fraction. Recoils can effectively de- 
crease the overall frequency of MBHs in small galaxies to 
~ 60%, while they have little effect on the frequency of 
MBHs in large galaxies (at most a 20% effect). 
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APPENDIX A: DEPENDENCE OF THE 
RESULTS ON THE MBH-HOST 
RELATIONSHIPS 

The estimates presented in the main text assume MBH 
mass es correlate with the velocity dispersion of their host 
halo l|Bandara et alj|2009h . and also that this relationship is 
redshift-independent . 

To test the robustness of our results against the as- 
sumed MBH-host relationship we implement a very different 
relationship, that is, we assume that the MBH mass scales 
with the mass of the dark matter halo: M = 10~ 6 Mh- This 
scaling leads to different dependencies of the MBH mass ra- 
tio, q, with the properties of the merging hosts (redshift, halo 
mass). The redshift dependence, however, cancels out in the 
expression for the MBH mass ratio, q, as long as the MBHs 
grow almost coevally. More important is the different scal- 
ing of q with the properties of the merging galaxies. In the 
"sigma-based" relationship we have: q oc (Mh,2/-Mh,i^)° , 
while for the "mass-based" relationship q oc (Mh,2/Mh,i). 
For Mh,2/Mh,i < 1 the sigma-based scaling leads to lower 
q values (this is true for any scaling relationship with an 
exponent larger than 1), and milder recoils. 

We find that the occupation fraction of small galaxies is 
mostly unaffected by the choice of the MBH-host relation- 
ship, and the difference for larger galaxies is within 20%, 
with the linear mass scaling leading to a lower occupation 
fraction. The redshift dependence of the MBH frequency is 
almost identical. 
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